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O ■ Abstract 
°. 

Using Woods-Saxon potentials and the eigenphase formalism for one-particle resonances, one- 

bJQ 1 

particle bound and resonant levels for neutrons as a function of quadrupole deformation are pre- 
sented, which are supposed to be useful for the interpretation of spectroscopic properties of some 

O , 

light neutron-rich nuclei with weakly-bound neutrons. Compared with Nilsson diagrams in text 
books which are constructed using modified oscillator potentials, we point out a systematic change 
of the shell structure in connection with both weakly-bound and resonant one-particle levels re- 
lated to small orbital angular momenta I. Then, it is seen that weakly-bound neutrons in nuclei 
such as 15_19 C and 33 ~ 37 Mg may prefer to being deformed as a result of Jahn- Teller effect, due to 
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' the near degeneracy of the ld 5 / 2 -2s 1 / 2 levels and the lf7/ 2 -2p 3 / 2 levels in the spherical potential, 

respectively. Furthermore, the absence of some one-particle resonant levels compared with the 
Nilsson diagrams in text books is illustrated. 
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I. INTRODUCTION 



The study of one-particle motion in spheroidal potentials, which is the basis for the 
understanding of deformed nuclei, started in the fifties [l[ \l , [], |j] . In particular, the work by 
S. G. Nilsson [3] played an important role for years in providing the basis for the classification 
of experimental data on the spectra of stable odd-A deformed nuclei |5|. Since the nucleon 
separation energy in stable nuclei is 7-10 MeV, the spectroscopic analysis around the ground 
state of those nuclei has been successfully performed in terms of harmonic-oscillator wave- 
functions. In contrast, the recent study of nuclear structure close to the neutron drip line 
points out the unique role of weakly-bound neutrons with small angular momentum I and 
the importance of coupling to the nearby continuum of unbound states. Due to the absence 
of the Coulomb barrier, weakly-bound neutrons with small I have an appreciable probability 
to be outside the core nucleus. Thus, those neutrons are insensitive to the strength of the 
potential provided by the well-bound nucleons in the system. In particular, the behavior of 
S1/2 neutrons is an extreme example due to the absence of the centrifugal barrier for the £=0 
orbit, compared with weakly-bound large i neutrons, of which wave functions stay mostly 
inside the potential. The i{l + 1) dependence of the height of the centrifugal barrier affects 
sharply the presence (or absence) of one-particle resonant levels at a given positive energy. 
Taking spheroidal Woods-Saxon potentials, the properties of both weakly-bound neutron 
orbits and one- neutron resonant levels are studied in Refs. j§[ Q, lsj]. 

There have been a number of self-consistent Hartree-Fock (HF) calculations of light 
neutron-rich nuclei, in many of which deformed HF solutions are indeed obtained. (As 
a classic work we may mention Ref. ij.) However, almost all those calculations are done 
either using the expansion in terms of harmonic oscillator bases or confining the system 
in a finite box. To our knowledge, no deformed HF calculation is yet available, which is 
carried out by integrating in mesh of space coordinate with proper asymptotic behavior for 
r = Rmaxi at which the nuclear potential is negligible. If deformed HF calculations are 
carried out in the latter way, the numerical results are totally independent of the values of 
R m ax and, furthermore, it becomes possible to estimate one-particle resonant levels without 
any ambiguity. 

The effective interactions to be used in HF calculations of nuclei far away from the 
stability line are not yet properly fixed. Thus, if the interaction chosen is not appropriate 
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for those nuclei away from the stability line, the answer is not reliable. Moreover, one 
HF solution gives one self-consistent deformation, which is determined by the two-body 
interaction selected, assuming that all technical problems in obtaining HF solutions are 
solved. Thus, it is not easy either to pin down the origin of the deformation obtained or to 
evaluate the ambiguity coming from the choice of the two-body interaction. 

In the analysis of observed spectroscopic properties of light neutron-rich nuclei with 
weakly-bound neutrons the shell model is so far used in most cases. The shell model should 
be applicable to those nuclei, if the configuration space is sufficiently large and the weakly- 
bound particles are properly treated. The latter condition is usually not satisfied, since 
harmonic-oscillator wave-functions in a limited space are used in most cases. The work 
summarized in Ref. [h]] is an exception among the calculations that can be systematically 
compared with experiments, however, the complicated calculations have so far been carried 
out only for the helium and oxygen isotopes. In any case, some physically interesting 
quantities such as one-particle energies (or shell structure) or nuclear shape (or deformation) 
are not directly obtained from shell-model calculations. 

Recognizing how useful it was to have Nilsson diagrams in the study of stable deformed 
nuclei |5( , which are constructed using modified oscillator potentials, in this paper we present 
some "Nilsson diagrams" that are relevant especially to some light neutron- rich nuclei to- 
wards the neutron-drip-line. Taking the Woods-Saxon potentials with the parameters ad- 
justed to some particular nuclei, both bound and resonant one-neutron levels are calculated 
as a function of quadrupole deformation. The change of nuclear shell-structure for neutrons 
is seen in both negative and positive one-particle energies of the Nilsson diagrams. The 
change comes from the unique behavior of neutron orbits with small i values, in particular 
£=0 and 1. The modified shell-structure has direct relevance to the ground and low-lying 
states of neutron-drip-line nuclei, in which weakly-bound neutrons are present. Considering 
the possible absence of many-body pair-field in light nuclei, the study of the present type 
of Nilsson diagrams can definitely help us to understand the origin of possible deformation 
and the related spectroscopic properties of light neutron-drip-line nuclei. 

In Sect. II some points of our model are summarized, while numerical results are presented 
in Sect. III. Conclusions and discussions are given in Sect. IV. 
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II. MODEL 



The occupancy of weakly-bound one-particle levels has a contribution especially to the 
tail of the self-consistent potentials. However, even for light nuclei presently considered the 
number of weakly-bound neutron(s) is much smaller than that of well-bound core nucleons. 
Namely, the major part of the nuclear potential is provided by well-bound nucleons. Thus, 
'or simplicity, the parameters of Woods-Saxon potentials are taken from the standard ones 



ll| for stable nuclei except for the depth, Vws- Namely, the diffuseness, the strength 



of spin-orbit potentials and the radius parameter are taken from those on p. 239 of Ref. 



The depth is adjusted so that a particular one-neutron level obtains a given desirable 



binding-energy in respective examples. 

The way in which bound one-particle levels are calculated is described in Ref. {(J, while 
the eigenphase formalism that is used to estimate one-particle resonant levels for a deformed 
potential is given in Refs. [?], Q] ■ The essential point is : the coupled equations obtained from 
the Schrodinger equation are solved in coordinate space with the correct asymptotic behavior 
of wave functions for r —>■ oo. The solution obtained in this way is totally independent of 
the upper limit of radial integration, R maxi if both the potential and the coupling term are 
already negligible at r = R ma x- One-particle resonant energy for [3 ^ is defined as the 
energy, at which one of eigenphases increases through it/2 as energy increases. In the limit 
of P — > this definition in the eigenphase formalism is in agreement with the definition of 
one-particle resonance in spherical potentials described in text books 12|; the phase shift 
increases through 7r/2 as energy increases. 

One-particle resonance is not obtained if none of calculated eigenphases do not increase 
through 7r/2 as energy increases. For example, we have no corresponding resonance in the 
case that a calculated eigenphase starts to decrease before reaching 7r/2 as energy increases. 
Even if one fails to obtain one-particle resonance defined in terms of eigenphase, for a 
certain small region of energy just after the disappearance of resonance the concentration of 
the wave functions inside the potential may still be found. However, the concentration will 
easily disappear after a short time if a resonance is no longer obtained in the eigenphase 
formalism. This situation is analogous to the case of the spherical potential, in which the 
phase shift starts to decrease before reaching ir/2 as energy increases [l^ . 

Compared with the Nilsson diagram based on modified oscillator potentials, the striking 
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difference of the level scheme obtained in the present work comes from the behavior of levels 
with low £ values (in particular, £=0 and 1) for (3=0 and those with small Q values (mainly 
Q n = l/2 + , l/2~ and 3/2~) for (3 ^ 0, in both the weakly-bound and positive-energy 
region. Note that the minimum £ value of possible components of Q n = l/2 + , l/2~ and 
3/2~ levels is equal to 0, 1 and 1, respectively. The absence of the centrifugal barrier for the 
£=0 channel produces the unique behavior of weakly-bound and positive-energy Q 71 = l/2 + 
orbits. However, we find that some Q n = l/2 + resonant levels survive in a higher-energy 
region (see, for example, the [200 1/2] level in Fig. 1), if the relative probability of S1/2 
component inside the potential is smaller than a certain critical value jsj]. Since the height 
of the centrifugal barrier becomes lower for a larger nuclear radius, the unique behavior of 
£=1 components will be more easily seen in nuclei with larger mass. 

III. NUMERICAL RESULTS 
A. Neutron-rich C-isotopes 

Taking Vws — —40.0 MeV and the radius parameter for A=17, at (3 = in Fig. 1 we 
obtain £(lc? 5 / 2 ) = —560 keV and £(2si/ 2 ) = —415 keV. The resonant energy of the ld 3 / 2 
level is 5.60 MeV. The near degeneracy of the lds/2 and 2si/2 levels compared to the high- 
lying ld3/2 level exhibits that for the spherical shape the neutron number N=16 may behave 
like a magic number. For (3 ^ the Si/ 2 , d 3 / 2 , d 5 / 2 , g7/2 and g 9 / 2 channels are included 
in the calculation of positive-parity levels, while the px/ 2 , p3/ 2 , fs/2 and f 7 / 2 channels for 
negative-parity levels. The parameters are chosen approximately for neutrons in the nucleus 
g 7 Cn, since the observed neutron separation energy of 17 C is —730 keV and the nucleus 
is presumably prolately deformed. Examining the Nilsson diagram in Fig. 1, it is easily 
imagined that a few neutrons occupying the weakly-bound almost-degenerate lds/ 2 - 2si/2 
shells at (3 = may prefer to being deformed in order to gain the total energy (Jahn- Teller 
effect). This may be the case of neutrons in 15 > 17 ' 19 C ; of which the neutron separation energy 
is small. It is also noticed that the observed ground states of J 7 Cn and g 9 Ci3 with J 77 = 
3/2 + and l/2 + , respectively, may be in a natural way interpreted as the band-heads of 
the intrinsic [211 3/2] and [211 1/2] configurations for prolate deformation, when the level 
scheme in the Nilsson diagram is applied to those nuclei. For some experimental evidence 
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for the deformation of those C-isotopes, see for example [13I, Q • 

The [220 1/2] resonant level is not obtained for (3 < —0.12 and Eq > 0.22 MeV, because 
the predominant component of the [220 1/2] wave function inside the nuclear radius is S1/2 
and thus decays very quickly Due to the same reason, the continuation of the [200 1/2] 
resonant level cannot be found for (3 > 0.46 and < < 1.82 MeV, while for smaller (3 
values the predominant component of the [200 1/2] level inside the nuclear radius is d 3 / 2 
and thus the resonant level is sufficiently well defined. 

In order to illustrate the near degeneracy of the 2si/ 2 and ld 5 / 2 levels at (3=0 in Fig. 1, 
in Fig. 2 the energy eigenvalues of Woods-Saxon potentials are shown, which are obtained 
by varying the depth while keeping other parameters the same as in Fig. 1. The 2si/ 2 level 
crosses with the lds/ 2 level for | Vws I slightly smaller than 40 MeV. For reference, the value 



of V ws that is obtained by applying Eq. (2-182) of Ref. [llj to fC u is -41.3 MeV. 

If a larger diffuseness, a > 0.72 fm, is used in Fig. 1 keeping other parameters of the 
Woods-Saxon potential unchanged, for (3=0 the 2si/ 2 level appears lower than the ld.5/2 
level. Nevertheless, in the region of an appreciable size of deformation the structure of the 
Nilsson diagram coming from the 2si/ 2 and lds/ 2 levels remains approximately the same. 

B. Neutron-rich Mg-isotopes 

Taking Vws — —40.0 MeV, in Fig. 3 the Nilsson diagram for neutrons is plotted for 
the radius appropriate for A=31. For the parameters of Fig. 3 the 2si/ 2 level for (3=0 is 
well-bound and, therefore, it lies approximately in the middle of the 2d 5 /o and 2d 3 / 2 levels, 
just as obtained from the level scheme in the modified oscillator potential [5j , but in contrast 
to the near degeneracy of the 2si/ 2 -ld 5 / 2 levels shown in Fig. 1. For (3 7^ the Si/ 2 , d 3 / 2 , 
d 5 / 2 , g7/ 2 and g 9 / 2 channels are included in the calculation of positive-parity levels, while the 
P1/2, P3/2, f5/2, h/2, hg/ 2 and hn/ 2 channels for negative-parity levels. The value of Vws is 
chosen so that the spin-parities of the ground state of nuclei ^Mg^ and ^ 2 Mg 21 , l/2 + [15] 
and 3/2 



16j, respectively, are obtained for (3 ~ 0.5 h may be identified 

as the band-heads of the [200 1/2] and [321 3/2] configurations with observed energies of 
about -2 MeV. 

The U7/2 resonant level at (3=0 is obtained at 6^=1.59 MeV. In contrast, the 2p 3 / 2 and 
2pi/ 2 resonant levels are not obtained, while the trace of the 2p 3 / 2 resonant level is expected 
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to lie below the lf 7 / 2 resonant level (see Fig. 4). Nevertheless, the energy is too high for 
the existence of a resonance with 1=1. The complicated behavior of the energies of the 
= l/2~ and 3/2~ levels for (3 < and below 2 MeV in Fig. 3 indicates the influence of the 
2p 3 / 2 and 2pi/ 2 levels which do not appear as resonant states around (3=0. Or, more exactly 
speaking, taking the fl n = l/2~ level expressed by the dotted curve as an example, the slope 
in the region of (3 < —0.3 indicates the 2p 3 / 2 level at (3=0 lying lower than the lf 7 / 2 level, 
while the slope for —0.2 < (3 < —0.14 suggests the 2px/ 2 level at (3=0 lying higher than the 
If 7/2 level. 

In order to locate the trace of the 2p 3 / 2 and 2pi/ 2 levels for Vws = —40 MeV in Fig. 
3, in Fig. 4 the energy eigenvalues of Woods-Saxon potentials are shown as a function of 
Vws while keeping other parameters the same as in Fig. 3. The 2p 3 / 2 and 2p!/ 2 resonant 
levels are not obtained for e > 0.80 MeV and e > 0.62 MeV, respectively, while the £=3 
one-particle resonant levels continue to be well defined above e=2 MeV. From Fig. 4 it is 
seen that the 2p 3 / 2 resonant level crosses with the lf 7 / 2 resonant level around V ws = —42 
MeV. Since in the present parameterization of Woods-Saxon potential the strength of the 
spin-orbit potential is proportional to Vws, it is somewhat misleading to draw the figure like 
Fig. 4 for a large variation of Vws- Nevertheless, using Fig. 4 it is easy to locate the trace 
of the 2p 3 / 2 resonant level for Vws = — 40 MeV, which is the depth used in Fig. 3. 

The fact that the trace of the 2p 3 / 2 resonant level is expected to lie below the lf 7 / 2 level 
in the positive energy region of Fig. 3 indicates that N=28 is not a magic number in the 
example. The near degeneracy of the lf7/ 2 and 2p 3 / 2 levels in the positive-energy region 
is indeed similar to the level scheme of almost degenerate ld 5 / 2 and 2si/ 2 bound levels in 
Fig. 1. This near degeneracy gives certainly the origin of possible deformation when a few 
weakly-bound neutrons occupy the lf 7 / 2 - 2p 3 / 2 shell. Namely, the degeneracy can be used 
to take a particular combination of the components so as to lower some level energy as 
deformation sets in (Jahn- Teller effect). This situation may correspond to neutron-rich Mg 
isotopes. 

In Fig. 3 the level coming from the 2p!/ 2 level is totally missing since it does not exist 
as a resonant level. The level denoted as [321 1/2] is obtained as a resonant level only for 
(3 > 0.24 and < 1.67 MeV. For (3 < 0.24 the energy of the level becomes larger than 1.67 
MeV, and the width becomes extremely large because £=1 is the predominant component 
of the wave function inside the nuclear radius. Therefore, the level cannot be identified as a 
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resonance. The one-particle resonant level obtained for (3—0, which lies next lowest to U7/2, 
is the lf 5 / 2 level found at 8.96 MeV that lies outside the range of Fig. 3. 

Taking V ws = —40.0 MeV, in Fig. 5 the Nilsson diagram is plotted for the radius 
appropriate for A=37. The H5/2 and 2p 3 / 2 resonances at (3=0 are found at 5.22 and 0.018 
MeV with the widths 2.08 and 0.005 MeV, respectively. Since e(lf 7 / 2 ) = -0.66 MeV, the 
distance between the lf 7 / 2 and 2p 3 / 2 levels is 680 keV, which is again very small compared 
with the distance obtained in the case that both levels are well-bound. This near degeneracy 
of the lf 7 / 2 and 2p 3 / 2 levels at (3=0 suggests that weakly-bound neutrons in nuclei with 
N=21-26 may prefer to being deformed. In the calculation of Fig. 5 for (3 7^ the si/ 2 , d 3 / 2 , 
^5/2) g7/2 and gg/ 2 channels are included in the calculation of positive-parity levels, while 
the pi/2, p 3 / 2 , h/2, $7/2, hg/2 and h n / 2 channels for negative-parity levels. The value of Vws 
is chosen so as to simulate the possible neutron-drip-line nucleus f 2 Mg 25 , which may have a 
neutron separation energy of a few hundreds keV. 

The 2pi/ 2 resonant level is not obtained at (3=0, while for (3 7^ no £l n = 1/2" one- 
particle level connected to the possible 2pi/ 2 level can survive as a resonant level. The fi 71 " = 
3/2~ resonant level (denoted by the dashed curve in Fig. 5) connected to the 2p 3 / 2 level at 
(3=0 cannot survive for (3 > 0.21 and En > 1.31 MeV, because the predominant component 
of the wave function inside the nuclear radius is £=1 and the level decays out quickly due 
to the low centrifugal barrier. The £l n = l/2~ resonant level (denoted by the dotted curve 
in Fig. 5) coming from the lf 5 / 2 level cannot survive as a resonance for (3 > 0.12 due to the 
increasing £=1 component inside the nuclear radius. 

In the region of a few MeV excitation-energy of Fig. 5, for both spherical and prolate 
shape we find no well-defined one-neutron resonant levels except the [303 7/2] level, since 
Nilsson levels expected in the region have either Q 71 = 1/2" or 3/2". 

It is noted that the energy distance between the lf 7 / 2 and 2p 3 / 2 levels at (3=0 in the 
Woods-Saxon potential becomes as large as several MeV when both levels are well bound, 
as known from the presence of the magic number N=28 in stable nuclei. 

IV. CONCLUSIONS AND DISCUSSIONS 

A few examples of Nilsson diagrams with both bound and resonant levels are given, the 
parameters of which are chosen to be appropriate for some light neutron-rich nuclei with 



8 



weakly-bound neutrons, using Woods-Saxon potentials. The absence of centrifugal barrier 
(very low centrifugal barrier) for £=0 (£=l) neutrons produces the 2si/ 2 (2P3/2) level close 
to the lds/2 (lfy/2) level, for both weakly-bound and low- lying resonant neutrons. This near 
degeneracy of the 2si/ 2 -ld 5/ / 2 and 2p3/ 2 -lf 7/ / 2 levels at /5=0 is recognized as the basic element 
of producing deformation for some neutron-rich C-Mg isotopes, in the case that the proton 
configuration allows the deformation. 

One-neutron resonant levels for (3 7^ are estimated using the eigenphase formalism. The 
Q w = l/2 + resonant level can hardly survive when the predominant component of the wave 
function inside the potential is si/ 2 , while the Q w = l/2~ and 3/2~ resonant levels are not 
obtained if the predominant component has £=1 and the energy is higher than 2 MeV in 
nuclei with A>16. In some nuclei the absence of those Q n = l/2 + , 1/2 - and 3/2~ resonant 
levels produces a low density of one-neutron resonant levels in the region of several MeV. 
How much the low density affects the many-body correlation such as pair correlation in 
nuclei towards the neutron drip line is a future problem to be studied and may be properly 
studied, only when the many-body correlation is studied treating the nearby continuum in 
a reasonable manner without discretizing the spectra. It is noted that the neutron one- 
particle levels obtained from Nilsson diagrams for (3 7^0 are those counted for band-head 
states of odd-N nuclei. Thus, rotational states, which are constructed based on those band- 
head states, should be in principle observed using a proper experimental method though 
high-spin states will have very narrow widths in the low-energy region. 

If one fails to treat properly weakly-bound neutrons or low-energy neutron resonant levels 
with small £, the HF 2p 3 / 2 (2si/ 2 ) level will not come down close to the lf 7 / 2 (ld 5 / 2 ) level. In 
any case, if the lf7/ 2 or ld 5 / 2 level is appreciably isolated, in the absence of pair correlation 
it may be possible to obtain an oblate shape as the deformation of the system with a few 
neutrons in the lf7/ 2 or lds/ 2 shell. This is because a preferred deformation is oblate at the 
beginning of the shell filling if a single j-shell is isolated (for example, see 20|, |21[), while it 
is prolate if shells with different j values are nearly degenerate as in the harmonic-oscillator 



potential 



221 ] . It is an interesting question whether any oblate shape is observed around the 



ground state of light neutron-drip-line nuclei of C-Mg isotopes. 
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Figure captions 



Figure 1 : Neutron one-particle levels as a function of quadrupole deformation. Parameters 
of the Woods-Saxon potential are designed approximately for the nucleus 17 C. One- 
particle levels are denoted by the asymptotic quantum numbers [N n 2 A Q]. The 
Q- values are denoted for four positive-parity levels for (3 < 0, since it may be difficult 
to see the connection to the levels for (3 > 0. One-particle levels appearing at (3=0 are 
lpi/2, ld 5 / 2 , 2si/ 2 and ld 3 / 2 levels at —6.77, —0.56, —0.42 and +5.60 MeV, respectively. 
One-particle levels in the positive-energy region, of which the phase shift (one of 
eigenphases) for (3 = {(3 ^ 0) does not increase through 7r/2 as energy increases, 
are not plotted. The neutron numbers 8 and 16, which are obtained by filling in all 
lower-lying levels, are indicated with circles. 

Figure 2 : Neutron one-particle levels as a function of the depth of Woods-Saxon potential, 
V W Si for (3=0. All parameters other than V ws are the same as in Fig. 1. Note that 
Vws — ~ 40 MeV is used in Fig. 1. The £=2 one-particle resonant level continues to 
be well defined above e = 2 MeV. 

Figure 3 : Neutron one-particle levels as a function of quadrupole deformation. Parameters 
of the Woods-Saxon potential are designed approximately for the nucleus 31 Mg. The 
Q 77 = l/2~ levels are denoted by dotted curves, the 3/2~ levels by dashed curves, the 
5/2~ levels by dot-dashed curves and the 7/2~ levels by dot-dot-dashed curves, while 
positive-parity levels are plotted by solid curves. One-particle levels appearing at (3=0 
are ld 5 / 2 , 2si/ 2 , ld 3 / 2 and lf 7 / 2 levels, from the bottom to the top. Neither 2p 3 / 2 nor 
2pi/2 levels are obtained at (3 = as one-particle resonant levels and, thus, they are 
not plotted in the figure. The next low-lying one-particle resonant level for (3=0 is the 
lf 5 / 2 level at 8.96 MeV that lies outside the range of the figure. See the text for details 
and the caption to Fig. 1. 

Figure 4 : Neutron one-particle levels as a function of the depth of Woods-Saxon potential, 
Vws, f° r (3=0. All parameters other than Vws are the same as in Fig. 3. Note that 
Vws = — 40 MeV is used in Fig. 3. The £=3 one-particle resonant levels continue to 
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be well defined above e = 2 MeV, while the 2p 3 / 2 and 2p^ 2 resonant levels do not 
survive for e > 0.80 MeV and e > 0.62 MeV, respectively. 

Figure 5 : Neutron one-particle levels as a function of quadrupole deformation. Parameters 
of the Woods-Saxon potential are designed approximately for the nucleus 37 Mg. One- 
particle levels appearing at (3=0 are 2si/ 2 , ld 3 / 2 , lfr/2, 2p 3 / 2 and lf 5 / 2 levels at —7.02, 
-5.28, -0.66, +0.018 and +5.22 MeV, respectively. The 2p 1/2 level is not obtained as 
one-particle resonant level. See the caption to Fig. 3 and the text for details. 
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Neutron one-particle levels in Woods-Saxon potential 

R = 3.990 fm (A = 31 ) a = 0.67 fm 
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Neutron one-particle levels in Woods-Saxon potential 
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